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Heat kernel coefficients encode the short distance behavior of propagators in the presence of 
background fields, and are thus useful in quantum field theory. We present a Mathematica program 
for computing these coefficients and their derivatives, based on an algorithm by Avramidi [n| . 
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1. PROGRAM SUMMARY 



Title of program: HeatK 
Version number: 1.0 



Available at: http : //f ermi . pha ■ jhu . edu/~booth/HeatK/ 



Programming Language: Mathematica 2.2 or higher 
Platform: Any platform supporting Mathematica 2.2 
Keywords: Heat kernel, HDMS coefficients, Effective Action 
Nature of physical problem: Computation of heat kernel coefficients. 
Method of solution: The non-recursive solution given by Avramidi |jj 
Limitations: Degree of the heat kernel coefficient. 

2. INTRODUCTION 

Heat Kernel (also know as HDMS or Seeley-de Witt) coefficients have many applications in physics and mathematics 
( a comprehensive list can be found in a paper by Avramidi and Schwimming [El). Because of this interest, there is an 
extensive literature describing their computation and properties. Unfortunately for those who are primarily interested 
in using the coefficients as inputs to a calculation, there is not one comprehensive source for them, and moreover a 
number of different sign and normalization conventions are in use, making it difficult to combine sources. A further 
complication is that many applications require derivatives of the coefficients, of which only a few can be found in the 
literature. HeatK was created to resolve these problems and provide an environment which allows easy manipulation 
of the results. 

3. THEORY 

Consider a general manifold of dimension d with metric c^ t „, fields ip^ and derivative acting such that [V M , V 
= F^ip and [V M , V„]ip a = R a Tflv ip T . For a differential operator H = — □ +Q, with □ = j^VV and Q a matrix- 
valued potential, the heat kernel expansion is developed by considering the exponential U(t) — exp(— tH), which 
satisfies the equation 

^ - -num. m 

The kernel of the operator U is then parameterized asQ 

- T^AV^^i-^, A (2) 

v ' fc>0 

where aj, are the heat kernel coefficients, tr(a;, x') is half the square of the geodesic distance between the points x and 
x' and 

A(x, x') = g- 1/2 {x) det (-V^V^z, x')) g~ 1/2 {x'). (3) 
is the Van Fleck-Morette determinant. 



*We follow the conventions of Ref. |l| throughout. 
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To obtain expressions for the a k , we make use of the non-recursive algorithm of Avramidi . Here we content 
ourselves with only a brief outline of the method, and refer the reader to the original papers for a more detailed 
description. 

Given the above data, we can then construct the two tangent vectors, a v = V u a and a u > = ^7 u >cr, noting that 
\a v a v = a = \o V 'O v . We also have the coincidence limits 

M = K1 = K'] = 0, (4) 

with 

[f(x,x')] = lim f(x,x'), 

X — >x 

It is also natural to define the tensor -rf- „ = V^er^ 1 and its inverse 7 / V ■ We can then define the parallel displacement 
operator V(x,x') as the solution to a^V^V = with [V] = 1 (When acting on vector indices, it is conventional 
to write g^ u / instead of 7-"V)- ^ n what follows it is convenient to work with quantities bases at x' , so we define 
V V ' = 7W M and the "gauge" field = P^V^V. 
Introduce a functional basis 

|0) = 1, (5) 
\n) = \u l l ...u' n )= K —J-o<-.-o<, (n>l), (6) 



(m| ee fa • ■ -n' m \ = (-l) m <i ' -^ V (mi • •• V^KO, (7) 



(m|n)= / d n x{n' 1 ---n' m \v[...v' n ), (8) 



its dual, 



and the scalar product 



so that 

(m\n) = s mn = S mn ■•■S^ly (9) 

For a field tp we then have 

(m\<p) = [V^-'-V^jp] , (10) 

and the covariant Taylor series is simply 

\<p)=vj2\ n )( n \<p)- (ii) 

n>0 

The advantage of this machinery is that it allows us to take the familiar HDMS recursion relation: 

1 + 7<r"V M J a k = A- 1 / 2 i/(A 1 / 2 a fc -i) ee Fa k -i, (12) 



k" '" y 

which follows from inserting the expansion Eq. (|^) into Eq. ([!]) and give meaning to its formal solution 



1„\-^/ 1 



a k =V\l + -Dj F^l + ^-jDj F---(1 + D)- L F, (13) 
(where D = a • V and F = T'^ 1 FV). This is accomplished by noting that D\n) = n\n), so that 

n>0 



from which it follows that 
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with 



(n\a k ) = 



k 



k - 1 



k + n k — 1 + Uk-i 1 + ni 



(m\F\n) = 



ni,---,rifc_i>0 

x(n|^|nfc_i>(nfc_i|F|n fc _2) ■ ■ ■ (ni|F|0), 



-f-D™ / 



77. 



Tl — 1 y (Ml— Mn-1 Mn"'Mm) 



(15) 



where 



m 



c(lM"-I'n-2 yV„-lVn) 



77,-2/ (Ml'"Mn-2 Cn-l-Mm) 



= (- 1 )" 



A Mi-Mr, 



(16) 



= (-ir 

Here X, Y and Z are simply the coefficients of the Hamiltonian operator transported to the point x', that is 

F = V^{A 1/2 □ A~ 1/2 + Q)V 



= X» V V^V^+F" V„/ + Z. 



(17) 



3.1. Derivatives 

Avramidi's formalism is naturally adapted to computing the symmetrized derivatives of the at, but frequently one 
requires the un-symmetrized derivatives. In order to obtain these, we can continue in the spirit of Avramidi and 
introduce matrix elements for the derivative operators. Defining 

TZi:fi m [h Xl ... Xn } = <0|V M1 . .. V„ m Vh Xl ... Xn \n), (18) 

we then have 

[V,, . . . V„ n a fe ] = £ Tfc-ft [l](m\a k ). (19) 

m 

Noting that 

V nVh Xl ... Xl \n) = -Vh Xl ... Xn r, { ; l \n- +V V v \(A v >h Xl ... Xn +^2& r Xi wh Xl ...a...x n )\n), (20) 

i 

where B is a special case of A: B a > u M = g° r a '^/ f ^g v a , it is possible to derive a recursive expression for T: 

m — 7i— 1 

T <-"' n r. i__ V- T K-<l«i-<lr i/J) ft i 

J Ml ■■■Mm [""Al •••A; J ^/^ -t(Ul.../i m _l L'/ Mm(Ki---«fe)'"'Al...A ! J 

fe=0 
m—n—1 

fc=l 
m—n—1 i 

+ E E T ^--^-V' K ' fc [^%™(«^--<)^---^]- ( 21 ) 



fe=l i=l 
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4. PROGRAM DESCRIPTION 
4.1. Functions 

The program implements most of the above functions, with the following correspondence: 



[V(, 



L J 


= DCLfkl, 

L J " 






= DCL[k, /Lti, . . 


• ,Mn], 


xf . 

(fl. ..!/„) 


= XX[a,/?,{"i 


. ■ • ■,"*}], 


vo 

y (^l...^n) 


= YY[a,{u u .. 


■,^n}], 




= ZZ[{l/!, 




•■V Mn) Q] 


= QQ[{A*i 


ML 


[Q] 


= QQ[ ], 






= f[m, H> 






= R[jti,i/,a,/3]. 





In general, as already seen above, we use a list of bracketed indices to indicate symmetrized differentiation, so we also 
have for example 

V (M1 . . . V M „)F^ - = A {Mi, • • • , M»}]. (23) 

For completeness, we also list the following functions, although they are essentially internal functions: 

(m\F\n) = MF[{z/ l5 . . . , u n ], . . . , A*m>], 

t;i:::Z [^...aJ = kd[{v u ■ ■ ■ , M. ■ ■ ■ . {Ai . . . A n }, h Al ... A j, (24) 

where we note that contrary to the convention, MD is not symmetric in the first set of indices. 



4.2. Options 



For the primary function DCL there are several options which control the normalization and degree of simplification. 
They are: 

HeatKNorm 

This controls the normalization that is used for the heat kernel coefficients returned by DCL. The default is 
Avramidi, which naturally invokes Avramidi's normalization (— l) k /k\. It can also be set to a function which 
should expect two arguments: the degree of the coefficient, and the number of derivatives. 

EvaluateK 

The curvature tensor occurs frequently through its symmetrized derivative K a p( n ) = V( Vl ...„ n _ 2 R a v n _ 1 \/3\ Vn )- 
Setting this option to True expands K in terms of R. 

EvaluateSymmetricDerivs 

Setting this to True means that objects such as Q[{^i, • ■ • , f n }] will be expanded in terms 

of explicitly symmetrized derivatives (represented by the operator Dx), so that for example, 
Q[{Mi,M2}] = (Dx[Q[ ],Mi,M2]+Dx[Q[ ], /i 2 , Ml ])/2. 

DoSimplif ications 

This option controls whether derivatives will be re-ordered and other simplifications performed. 
DoBasicSimplif ications 

This option controls whether to do basic simplifications like contracting indices. Setting this to False is really 
only useful when debugging. 

There is also a global option defined for the HeatK: 
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EvaluateZZ 

Setting this option to False means that ZZ, YY and XX will not be evaluated. This makes the results much more 
compact, though cryptic. 

Two variables which describe the form of the derivative operator. The first, $HeatKMetricStyle indicates whether 
one wishes to use the GR style metric diag(—l, 1, . . . , 1) or the HEP style metric diag(l, —1, . . . , —1). Setting it 
to —1 (the default) selects the former. The second variable, $ScalarCurvatureFactor is used to include an ex- 
plicit factor of the scalar curvature (R) into the Hamiltonian. For example, with $HeatKMetricStyle = — 1 and 
$ScalarCurvatureFactor = —1/6, one has H = — □ +Q — R/6. 

4.3. Examples 

The following presents a simple example illustrating the use of the package, and showing how to change the settings 
to conform to various conventions. First load the package, assuming it is in the search path: 

In[l] := «HeatK.m 

By default, the package is set up to follow the sign and normalization conventions of Avramidi's article Q We can 
check that it reproduces those results. 

In [2] := al = DCL[1] 

Out [2] = ZZ 

;() 

In [3] := a2 = DCL [2] // Expand 
Out [3] = ZZ 

; (ul.ul) 

ZZ ** ZZ 

;() ;() 3 

In [4] := a3 = DCL [3] // Expand 

0ut[4]=-YY ** ZZ YY ** ZZ 

ul;(u2) ;(ul,u2) ul;(u2,u2) ; (ul) 



ZZ ** ZZ ZZ ** ZZ 

;() ;(ul,ul) ;(ul) ; (ul) 



ZZ ** ZZ 

;(ul,ul) ;() 

+ ZZ ** ZZ ** ZZ + 

2 ;() ;() ;() 

XX ZZ ZZ 

ul,u2; (u3,u3) ;(ul,u2) ; (ul ,ul ,u2 ,u2) 

+ 

6 10 



It is not difficult to see that the above results agree with those presented by Avramidi. Now, if we prefer, we can 
follow the conventions in his thesis H. 

In [5] := $HeatKMetricStyle = -1; 
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In [6] := al = DCL[1] 
Out [6] = ZZ 

;() 

In [7] := a2 = DCL [2] // Expand 
Out [7]= ZZ 

; (ul.ul) 

ZZ ** ZZ + 

;() ;() 3 

In [8] := a3 = DCL [3] // Expand 

0ut[8]=-YY ** ZZ YY ** ZZ 

ul;(u2) ;(ul,u2) ul;(u2,u2) ; (ul) 
+ 

3 4 

ZZ ** ZZ ZZ ** ZZ 

;() ;(ul,ul) ;(ul) ; (ul) 
+ + 

2 2 

ZZ ** ZZ 

;(ul,ul) ;() 

+ ZZ ** ZZ ** ZZ + 

2 ;() ;() ;() 

XX ZZ ZZ 

ul,u2; (u3,u3) ;(ul,u2) ; (ul ,ul ,u2 ,u2) 

+ 

6 10 

Now let's now back up and see what some of these expressions are in more familiar notation. We can do this by 
changing the options for HeatK, allowing XX, YY and ZZ to be evaluated. We also switch back to the conventions in 
the nuclear physics article 0] . 

In [9] := SetOptions [HeatK , EvaluateZZ -> True] ; 
$HeatKMetricStyle = -1; 

In[10]:=ZZ[] // ContractAll 

Out [10]= R 

QQ — 

6 

Simplif yCurvature uses some tricks to put the curvature tensor into a "canonical" form, while ContractAll and 
ReduceContractions perform contractions, and simplify them - ZZ, XX and YY do not simplify their results , because 
they are usually not called directly. 

In[ll] := ZZ[ul, ul] // ContractAll // Simplif yCurvature // 
ReduceContractions 

Out [11]= -F ** F R 

ul,u2 ul,u2 ;(ul,ul) 

+ QQ + 

2 ;(ul,ul) 5 



G 



F R R R 

ul,u2 ul,u2 ul,u2 ul,u2,u3,u4 

+ 

3 30 30 

We see that these results agree with Eqns. 3.83a and 3.83b of Ref. [l| 

Now let's compare with a standard reference, Barvinskii and Vilkoviskii's (BV) Physics Reports Q. Avramidi 
defines his heat kernel coefficients with an explicit weight of 1/fc!, so we have to divide by l/k\ to get the more 
standard coefficients used by BV. Also, BV uses (is) rather than t as the exponential parameter introduces an extra 
factor of (— l) fe in their definition of the coefficients. All told, their conventions correspond to the following: 

In [12] := SetOptions[DCL, HeatKNorm -> ( (-1) " (#1) / (#1) ! &)] ; 

$HeatKMetricStyle = -1; $ScalarCurvatureFactor = 1/6; 

We must also note that their matrix P is opposite in sign to Avramidi's QQ, so we correct for that as well. 
In[13] := DCL[1] /. _QQ -> -QQ 
Out [13]= QQ 

In [14]:= (a2 = DCL[2] // Expand) /. _QQ -> -QQ 

0ut[14]=D D (QQ) D D (R) F ** F 

ul ul ul ul QQ ** QQ ul,u2 ul,u2 

+ + + 

6 180 2 12 

2 2 
F R R R 

ul,u2 ul,u2 ul,u2 ul,u2,u3,u4 

+ 

36 180 180 

In [15] := DCL[1, ul] /. _QQ -> -QQ 

0ut[15]=D (QQ) D (F ) 
ul u2 ul,u2 

+ 

2 6 

In [16] := DCL[1, ul, u2] /. _QQ -> -QQ 

0ut[16]=D D (QQ) D D (F ) D D (F ) 

ul u2 u2 u3 ul,u3 ul u3 u2,u3 

+ + 

3 12 12 

D D (R) D D (R ) QQ ** F 

ul u2 u3 u3 ul,u2 ul,u2 

+ + + 

180 60 6 

F ** QQ F ** F F ** F 

ul,u2 ul,u3 u2,u3 u2,u3 ul,u3 

+ + 

3 12 12 

R R F R 

ul,u3 u2,u3 u3,u4 ul,u3,u2,u4 
+ + 
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45 



18 



R 



, R 
u3,u4 ul,u3,u2,u4 



F 



R 

u3,u4 ul,u4,u2,u3 



+ 



+ 



180 



18 



R 



R 

u3,u4 ul,u4,u2,u3 



R 



R 

ul,u3,u4,u5 u2,u3,u4,u5 



+ 



180 



90 



More Examples (as Mathematica notebooks) can be found at the WWW-site listed in the program summary. 
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